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ABSTRACT 

Galactic stars belonging to the Of?p category are all strongly magnetic objects exhibiting rotationally modulated 
spectral and photometric changes on timescales of weeks to years. Five candidate Of?p stars in the Magellanic Clouds 
have been discovered, notably in the context of ongoing surveys of their massive star populations. Here we describe 
an investigation of their photometric behaviour, revealing significant variability in all studied objects on timescales of 
one week to more than four years, including clearly periodic variations for three of them. Their spectral characteristics 
along with these photometric changes provide further support for the hypothesis that these are strongly magnetized O 
stars, analogous to the Of?p stars in the Galaxy. 
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1. Introduction 


The spectral type Of?p, defined bv IWalbc^ (|l972li . may 
at first sight seem awkward with its question mark, but 
it is not a typo at all. The “f?p” was introduced to iden¬ 
tify stars showing strong emission lines of N ill A 4634-41 
and He IIA 4686 (hence the “f”), with a spectrum similar 
to those of supergiants, but showing other different charac¬ 
teristics including the unusual presence of C ill A 4650 lines 
in strong emission, leading to their peculiarity flag (“p”); 
a question mark was then added to express doubt that 
these stars are normal Of supergiants. Only three Galactic 
objects were classified in that cate gory in th e seventies: 
HD 108, HD 191612, and HD 148937 TWalbornl. [TttI . 

In the last decade, a growing interest in this class has 
appeared, as recurrent spectral variations were identified 
in their Balmer and He i lines. Derived periods am ount to 
^7 d for HD 148937 (iNaze et al.l. l2008a ; Dlazee^^ 


538 d for HD 191612 (IWafboru et afl 120041: 1^77^ 0"^! | Wade et al. |, | 2012aD . 


properties (same period, simultaneous extrema) as the 
changes in the optical spectra. 

After a strong m agnetic field was detected in HD 191612 
(jPonati et al.Ll200'^ . the magnetic oblique rotator scenario 
was proposed to explain t hese peculiar properties. The 
monitoring of I Wade et al.l (l201lD confirmed the proposed 
scenario for that star. In this model, the stellar winds of 
massive stars are confined in the equatorial regions by a 
strong dipolar magnetic field and, as the star rotates, these 
regions come in and out of direct view, generating the ob¬ 
served recurrent variations found at all wavelengths. In this 
context, the recurrence timescale is the rotational period, 
and long rotation periods may then be ex plained by mag¬ 
netic braking over the lifetime of the star (|ud-Doula et al.L 
l2009f) . This implied that all Of?p stars are strongly mag¬ 
netic, and HD 108 and HD 148937 w e re indeed soon found 
[2(M l2nTTl [ml iMartins et ahl 


TiTTTTT to be SO (Hubrig et al 


l2007t Howarth et ahl 2007 ), and ^55yrs for HD 108 


( Naze et al .1. 120011: 1 iNaze et al. I . I2006l1 . A strong X-ray emis- 


At the same time, new Of?p stars were identified in the 
Gala xy, with very similar characteris tics (IWalborn et a^ 


and CPD—28°2561 (iHubrig etAD 12011 . 20131: period of 


sion, brighter than that of norinal O stars, w as also found l2010l1 : NGG1624-2 ( period of 158 d, Wade et ahl . 12012 
for these stars ([Naze et al.l . l2004120071 [^08al1 . with a mod¬ 
ulation of the high- energy emission wi th the same period as 
the optical changes (iNaze et ^.l2010h. Finally, o ptical pho¬ 
tometric changes (|Naz412004 [Barannikov 2007f) as well as 
variations o f line profiles in the UV range ( Marcolino et al.l . 

I 2 OI 2 LI 2 OI 3 II were also reported, and they display the same 


* FNRS Research Associate 

** Operated by AURA, Inc., under NASA contract NAS5- 
26555. 


73 d, I Wade etAri 1 201511 . It must be stressed that, while 
about ten O stars are now known to be magnetic, the only 
actual class of magnetic massive stars is the Of?p category. 
Spotting the typical spectral properties of Of?p stars there¬ 
fore constitutes a powerful, albeit indirect, way to identify 
a magnetic 0-type star. However, some stars with spec¬ 
tral types 03.5-5.5 actually display strong Giii lines nat¬ 
urally (i.e. in the absence of a magnetic field), which led 
to the definition of a separate spectral type dubbed “Ofc” 
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Table 1. List of our targets, with their position and their star and field IDs in ASAS and OGLE databases. Shifts applied 
to OGLE-II and OGLE-IV datasets to bring them in line with OGLE-III data are also provided (in bold). 


Star 

RA 

DEC 

ASAS 

OGLE-II 

OGLE-III 

OGLE-IV 

SMC159-2 

00h49m58.72s 

-73“19’28.4” 


SMC SC5#95303 

(A — —0.042 mag) 

SMC100.8#52967 

SMC720.27#273 

(A = -0.057 mag) 

2dFS936 

00h53m29.95s 

-72°41’44.5” 

005329-7241.6, 005330-7241.6, 
005331-7241.8, 005330-7241.7, 
and 005329-7241.4 

SMC.SC6#237339 
(A = -0.013 mag) 

SMC101.2#21946 

SMC719.11#68 
(A — —0.032 mag) 

AV220 

00h59m09.97s 

-72°05’48.3” 

005910-7205.8 
and 005909-7205.8 


SMC108.3#6501 

SMC725.16#45 
(A — -^0.006 mag) 

BI57 

05h01m08.59s 

-68°11’45.1” 

050107-6811.6, 050109-6811.8, 
050111-6811.6, 050106-6811.6, 
and 050108-6811.8 


LMC125.2#9 

LMC510.32#25733 
(A — -1-0.040 mag) 
and LMC532.08#12629 
(A — 4-0.026 mag) 

LMC164-2 

05hl3m49.88s 

-69°23’21.7” 

051351-6923.1, 051346-6923.1, 
051349-6923.4, 051351-6923.1, 
051350-6923.4, 051350-6923.1 

LMC_SC9#216949 
(A — —0.003 mag) 

LMC111.2#10462 

LMC503.15#28081 
(A — 4-0.013 mag) 

Notes: 2dFS 936 is also known under the name 

[MD2001] Anon 1; ASAS dataset 051346-6923.1 was discarded as it contains only two flag ’A’ points 


and in one of them the star appeared much brighter than in all other data. 


(|Walborn et al.l . 120101: iMartins fc Hillieil l2012^ . The pres¬ 
ence of strong G ill lines is therefore not sufficient for clas¬ 
sifying a star as Of?p. Moreover, it is now known that the 
G III emission lines in some Of?p spectra disappear entirely 
at certain phases, while in others they remain strong or 
weak at all phases. Hence, the identification of an Of?p- 
type star is now rather performed through a large body of 
evidence: peculiar spectral characteristics (e.g. in the opti¬ 
cal: emission lines narrower than absorption lines, emissions 
or P Gygni features within Balmer and Hei lines,etc.) that 
are often periodically variable, small recurrent variations in 
brightness, abnormally bright X-ray emissions, and strong 
dipolar magnetic fields (see e.g. lNaze et al.l . l2008bll . 

Like magnetic fields, metallicity is a key ingredient for 
many aspects in massive stars’ lives (formation, evolution, 
winds). It is therefore of utmost interest to study whether 
magnetism also occurs, and with what characteristics, at 
other metallicities. The closest laboratories for such tests 
are the Magellanic Clouds (MCs), and the best way to find 
likely magnetic stars is to locate Of?p stars. Indeed, other 
kinds of magnetic stars, e.g. active M dwa rfs or Ap stars 
(jMaitzen et al.L l200ll iPaunzen et al.l . l20Illl , might be eas¬ 
ier to analyse, but the much higher intrinsic brightness of 
O stars makes them better targets in the MCs. In fact, pe¬ 
culiarities typical of Galactic Of?p stars ha ve already been 
detected for five stars in the MCs: AV22 0 (IWalborn et al 


200C; 


2001 


Massey fc Duffvl. 200111 . 2dFS936 (iMassev fc Duffy . 
Evans et al.l . l2004ll . BI 57 (1. Howarth, private com¬ 


munication; Walborn_etal.in_OTeparation), SMC159-2 and 
LMC164-2 (|Massev et al.l . l2014li. A sixth one, LMC N82 
(or Br ey 3a) was proposed bv iHevdari-Malay e^fc 
1992). but it was classified as Ofn bv IWalborn et al 


( 20031) . Nothing more is currently known about these five 


stars, as their great distance makes the investigations dif¬ 
ficult. This paper provides the first step forward, with the 
analysis of photometric data of these stars, with the aim 
of identifying a key characteristic of magnetic Of?p stars: 
periodic variability consistent with a long (weeks-years) ro¬ 
tational period. Section 2 presents the observations used to 
do so, while Section 3 presents our results, and Section 4 
concludes this paper. 


2. Data 

In Galactic Of?p stars, periods range from a week to 
decades. Trying to pinpoint the recurrence timescale from 


random spectral observations can thus be very time- 
consuming. It is even more the case when the faintness of 
the object imposes long exposures to get spectra of good 
quality, as for MG objects. However, Galactic Of?p stars 
do not only display spectral variations, they also show pho¬ 
tometric ch anges, of the o r der o f 0.05 mag fo r HD 108 and 
HD 1 91612 (iKoen fc Ever! . 120021: lNaz4 l2004 iBaranniko^ 
120071) . Brightness variations are much easier to detect for 
MG objects since there have been several photometric sur¬ 
veys of these regions over the years. 

We thus s earche d the archives of the OGLE- H 
(|Udalski et al.l . ll997fFl and ASAS (|Poimanskil Il997li^ 
projects. We found data (in the /-band for OGLE-H, Id- 
band for ASAS) for the five candidate Of?p stars in the 
MCs. From these datasets, only the best quality data 
were kept (box “good quality” selected for OGLE-H when 
downloading data - we note, however, that data with flag 
down to F remain; only data with flag A were kept for 
ASAS). For OGLE-H, two photometric datasets are avail¬ 
able: one called DIA wh i ch refers to Differen ce Image 
Analysis (lAlard fc Luntonl . Il998t IWozniak . 2000I). and oii e 
called PSF which refers to PSF fitting ( Szvmanskil . |200§) . 
The analyses were performed on both sets of magnitudes: 
they gave similar results. As DIA data are less noisy (and 
known to be more reliable, see ISzvmanskil 120051) . the num¬ 
bers presented below refer to them. 

Further OGLE data, from OGLE-III and OGLE-IV 
campaigns, were also ava ilable. The th ird phase of the 
OGLE project (OGLE-III, IUdalskiir2003[) was realized dur¬ 
ing the years 2001-2009 using an 8-chip CCD mosaic. The 
photometry was obtained using the DIA method. OGLE- 
IV is the currently ongoing, fourth phase of the project. It 
started in 2010, using a huge, 32-chip CCD mosaic and the 
DIA reduction scheme. Analysing the OGLE datasets we 
noted small residual magnitude shifts between the three 
campaigns. To correct for them, we computed, for each 
dataset, the mean of all data, the median of the 10% largest 
data (i.e. 90th percentile), and the median of the 10% small¬ 
est data (i.e. 10th percentile). OGLE-III data are consid¬ 
ered as the best calibrated amongst the OGLE collabora¬ 
tion, and were then used as reference: the differences be¬ 
tween the derived mean and median values of OGLE-H/IV 
and their equivalent in OGLE-III data yield shifts that were 
applied to the OGLE-H/IV before analysis. We note that 


^ http://ogle.astrouw.edu.pl/ 

^ http://www.astrouw.edu.pl/asas/ 


2 






















































































Naze et al.: Periods for extraglactic Of?p stars 


Table 2. Properties of the detected periods, 
Fourier method. 


from the 







14.26 

- 

ID 

To-2 450 000 

period (d) 

amplitude (mag) 


14.28 


SMC159-2 

3598.706±0.049 

14.914±0.004 

0.00719±0.00015 

W 

14.3 

- 

2dFS936t 

3993.3±1.7 

1370±30 

0.01504±0.00011 

o 

14.32 

- 

LMC164-2 

4027.338±0.015 

7.9606±0.0010 

0.00926±0.00011 


14.34 

_ 

BI57* 

5302.2±1.2 

400.0±3.5 

0.00550±0.00010 





To (hence = 0) corresponds to minimum brightness; * not strictly 
periodic (see text); ^there are additional periodicities in 2dFS936 (see 

text). 


using means or any of the medians yields the same shifts. 
For reference, they are shown in Table [T] along with the IDs 
of the fields and stars in the different programmes. 

Although it is currently off-lin e, we managed to get ac- 
cess to the EROS-2 database isee lTisserand et ^1120071 for 
details). Data are available for three of our stars (AV220, 
2dFS 936, and LMC164-2). They cover a shorter time inter¬ 
val than OGLE data, from HJD-2 450 300 to 2 452 700, but 
can be used as a confirmation tool. Other photometric sur¬ 
veys have no available data for these objects (e.g. Catalina 
sky survey) or have no public data. 

On each of these datasets, we first applied a 
test for constancy and then several period search algo¬ 
rithms. They fall in four categories: (1) the Eourier al¬ 
gorithm adapted to spar se/uneven datasets (|Heck et all . 
119851: iGo sset et akl. 12001 . a met hod rediscovered recently 
bv IZechmeister fc Kiirstei] 20091 - th ese papers also note 


that the method of Scargld (|1982ll . while popular, is 
not fully corre ct, statistically), (2) two different s tring 
length methods (|Lafler fc Kinmanl . ll965t Rensonl.ll978D. fS) 
three binned ana l yses o f variances ( Whittaker fc RobinsonI 
11944 l.lurkevichl 1 19711 which is i dentical, with no bin 
overlap, to th e “pdm” method of IStellingwer j [197^ and 
iGnypensI 1987 - which is i d entica l to the “AOV” method 
of ISchw arzenberg-Gzernv (llOSOTll, and (4) c onditional 
entropy (ICincotta et al.l . 19991 : ICincottal . Il999l see also 
iGrahain et al.1 l2013ll . Each of these methods has its ad¬ 
vantages and its drawbacks - the most reliable being the 
Fourier method while the fastest ones usually are analyses 
of variances - but a period clearly identified by all of them 
is certainly real. All uncertainties reported in this paper 
correspond to Icr errors. 

In view of the results (see next section), additional spec¬ 
tra are needed. Indeed, usually, only one spectrum for each 
object exists in the literature, but more than one spectrum 
is needed to examine variations. Spectral monitorings of 
AV220, BI57, and 2dFS936 are examined in a compan¬ 
ion paper (Walborn et ah, in preparation). For the remain¬ 
ing two objects, the published spectra of LMC164-2 and 
SMG159-2 were o btained in the frame work of a WR sur¬ 
vey in the MCs (iMassev et al.l . l2014l) . This campaign is 
not completed yet, which enabled us to get one additional 
spectrum of SMG159-2 in September 2014. As it is the sole 
new one, it will be presented here rather than by Walborn 
et al. (in preparation) who pres ent larger spec t ral se ts. As 
with the two previous spectra (IMassev et al.l . l2014ll . this 
new spectrum was obtained at Las Gampanas Observatory 
with the 6.5-m Magellan II (Cl a y) tel escope using the MagE 
spectrograph (iMarshall et al.l . 120081) and a 1” slit which 
provided a resolution ranging from 0.7 to 1.7 A going from 
the blue to the red end of the optical spectrum. The ex¬ 
posure time was 3 x600s. The data were processed with a 
combination of the IRAE mtools package, originally devel¬ 



2000 4000 6000 

HJD-2 450 000. 


Fig. 1. Photometric variations of 2dES 936, observed in 
the I-band (OGLE, top - OGLE-II data are shown by red 
crosses, OGLE-III by blue triangles, and OGLE-IV by black 
dots), non-standard B-band (EROS-2, middle), and V-band 
(ASAS, bottom). It must be noted that the top panel spans 
only 0.10 mag while the middle and bottom panels span 
much larger ranges (0.4mag and 1.2mag, respectively). 

oped by Jack Baldwin for the reduction of MI KE data, and 
standa rd IRAE echelle tasks as described bv IMassev et al.l 
(l2012l) . ThAr lamps were used to derive wavelength solu¬ 
tions, and flux standards were observed on each night in or¬ 
der to flux calibrate the data. The two spectra of SMC159-2 
were normalized considering clean continuum windows. 

3. Results 

Before looking at each object in turn, a note about the data 
quality should be made. We expect the photometric varia¬ 
tions to be very small (in the Galaxy, variations no larger 
than 0.05mag are seen, see previous section). OGLE data 
typically have errors of 0.003-0.004mag while the errors on 
AsAS and EROS-2 data are ten times larger (0.04-0.05 mag 
for ASAS and 0.015-0.15mag with a median of 0.03mag for 
EROS-2). In addition, the extraction apertures are different 
(from about 1.5” in OGLE to a radius of 15” in ASAS, with 
EROS-2 in-between), leading to different degrees of con¬ 
tamination. It is therefore unsurprising that the different 
datasets do not produce exactly the same results. A good 
example is provided with 2dES936: the coherent variations 
detected by OGLE have an amplitude of 0.015 mag and this 
is confirmed in EROS-2 data (see below), but this signal is 
buried in ASAS noise as can be seen from a direct com¬ 
parison of data taken by the two surveys at similar epochs 
(Fig. [I]). Therefore, we definitely trust more OGLE data 
than ASAS data. 


3.1. AV220 

For AV220, the photometric data show a dispersion that 
is significantly larger than the error bars (e.g. Fig. [5J left 
panel), and are thus found significantly variable in a test. 
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AV 220 




H,ID-2450000. Frequency 


Fourier Period Search 



0 0.005 0.01 0.015 0.02 


Frequency (time unit ') 


Fig. 2. Left: OGLE-III and IV lightcurve for AV220 (colour scheme as in Fig. [!}. Middle: Comparison of the different 
period search methods, with tickmarks at the top indicating periods of 3846 d, 1299 d, and 1100 d. The presence of a 
signal is marked by a maximum for Fourier and AOV methods, but by a minimum for conditional entropy, Renson and 
Lafler & Kinman methods. Right: Fourier periodogram for the raw (top) and prewhitened (middle, for the best-fit Fourier 
period P=1299d) data, along with the spectral window (bottom). 


BI 57 




HJD--2450000. Frequency 


Fourier Period Search 

0 0.005 0.01 0.015 0.02 



Fig. 3. Same as Fig. [5] for BI57 (a best-fitting double sinusoid, with the two periods of 787 d and 400 d is added in 
red in the left panel; these periods are indicated by tickmarks in the middle panel; prewhitening in the right panel was 
performed for P=400d). The potentially problematic OGLE-IV dataset (LMC532.08) is shown in green. 


However, period search algorithms do not find a clear and 
unequivocal signal (Fig. [21 middle panel). With the ASAS 
dataset, most methods display some low-amplitude peaks 
around 10-15 d, but periodograms look overall like white 
noise. For OGLE data, power is found only at the lowest fre¬ 
quencies, but different methods disagree on the value of the 
best period (1299d, 3846d?), and prewhitening with these 
frequencies does not yield a flat periodogram (Fig. |21 right 
panel). This may come from a changing variation timescale, 
with a variation apparently more rapid in the OGLE-IV 
dataset. Indeed, when looking at photometry over time 
(Fig. 121 left panel), brightness maxima are seen around 
HJD =2 452 900, 2 455 200, and 2 456 600 while minima are 
seen around HJD =2 452 000, 2 454400, and 2 456 000 (we 
note, however, that 2 455 200 and 2 456 600 are at the begin¬ 
ning and end, respectively, of the dataset). To assess this 
difference, we perform period searches on the OGLE-III 
and OGLE-IV datasets separately: periodograms based on 
OGLE-III data favour “periods” of 2500 d or 3700 d, while 
those calculated with OGLE-IV data yield shorter values 
of 1200-1400 d. The change with time is thus not an effect 


of the annual gaps, but a true variation. The EROS-2 data 
produce similar results and could not provide a stronger 
constraint on the timescales. A secure period has thus not 
yet been identified for AV220, and further monitoring is 
needed to better understand the photometric variations of 
this star. 


3.2. BI57 

Regarding BI57, it must be noted that, in the LMC532.08 
field, the star is located close to the outer edge of the mo¬ 
saic, where the template image is of worse quality and the 
photometry calibration less reliable owing to vignetting. 
However, discarding this dataset from the analysis does not 
change significantly the following results. BI57 displays ob¬ 
vious and significant variability of its photometry (Fig. [31 
left panel). In period searches, ASAS data yield a marginal 
detection at a Fourier period of P = 96.7 ± 0.3 d, while 
OGLE data clearly favour a period of 400±3d (Fourier) or 
its double (787±14d, entropy & variance analyses). While 
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BI57 

14.08 -m^^-1-^-1- 



J I I L 


Fig. 4. Photometric variations of BI57, phased using the 
best-ht Fourier ephemeris (period of 400d, see Tabled]). 
The thick red line corresponds to the best-fit Fourier sinu¬ 
soid for that period, while the symbols for data are as in 

Fig. 13 


the photometry appears rather well phased with the lat¬ 
ter periods (e.g. Fig. Sj), it must be noted that the dis¬ 
persion in the phased diagram appears higher than for the 
three following Of?p candidates, pointing to incompatibil¬ 
ities of the OGLE dataset with a strict periodicity. As 
further evidence, a model considering the sum of two si¬ 
nusoids with periods 400 d and 787 d does reproduce the 
average magnitudes of each observing run, but sometimes 
fails to reproduce the detailed trends (increase/decrease) 
of the observed variations (see such a model in red in the 
left panel of Fig. 13). Furthermore, prewhitening by these 
periods does not yield a flat periodogram as it should in 
the case of a strict periodicity (e.g. Fig. [3 right panel). 
Analysing the three OGLE datasets separately confirms the 
presence of differences: for example, the 787 d peak is higher 
in the Fourier periodogram of OGLE-III data while the pe- 
riodograms based on the two OGLE-IV datasets favour the 
400 d peak. Therefore, while it is clear that some changes 
occur with periods of 400 d and 787 d, such periodicities 
cannot explain the full range of variations in BI 57. 

3.3. 2dFS936 

The case of 2dFS936 is very different. The data appear 
significantly variable (the star was in fact r eported as an 
“irregular variable” bv lKourniotis et id]l2014f ) ; a clear long¬ 
term modulation is readily detected by eye in OGLE data 
(Fig.[T]). It is indeed detected independently by the different 
period search methods (Fig. [3 left panel), with a best-fit 
period of 1370±30d (Table[3- Prewhitening by this period 
yields a flat periodogram while phased photometry results 
in very coherent variations (Fig. [3 middle and right pan¬ 
els), securing that detection. In addition, the spectroscopic 
changes of 2dFS 936 seem to occur in phase with the same 


period (see Walborn et ah, in preparation): this 1370d pe¬ 
riod clearly appears as the rotation period of a typical Of?p 
star. 

It should be noted, however, that the maximum 
at HJD ~2 453 250 and the following minimum at 
HJD ^2 454000 occur when the star was fainter than for 
other maxima and minima. Examining further the peri¬ 
odogram, we found that a second set of peaks exists, along 
with their daily aliases: in the interval 0.-4.5d“^, their 
frequencies are 0.48525d“^ & 0.51246 d“^, 1.48676d“^, 
2.48952d-\ 3.48952d-i & 3.49224d-\ and 4.49228d-^ 
With the current dataset, it is difficult to assess which 
frequency is the correct one: data points prewhitened for 
the 1370 d variation yield nicely phased variations for all 
frequencies larger than ld“^ (e.g. Fig. Ej). The origin of 
this additional variation may not be linked to the Of?p 
phenomenon, it may rather be related to pulsations or to 
eclipses in a binary system. In the latter case, the orbital 
period would then be the double of one of the periods de¬ 
tected by period searches and the lightcurve with two sim¬ 
ilar eclipses would indicate similar stars. However, such a 
period (0.4-4 d) would be very short for such massive ob¬ 
jects and a line-of-sight coincidence between 2dFS 936 and 
an eclipsing binary could then be required. Identifying the 
correct frequency and understanding its nature implies re¬ 
observing 2dFS 936, this time over several nights with a 
high sampling frequency to avoid any ambiguities. 

Formally, using a test, the EROS-2 data in their blue 
band (where most data points are available) do not deviate 
significantly from a constant, but period search methods, 
being very sensitive to periodic signals, readily confirm the 
signals found in OGLE data: main period at 1389±80d with 
an amplitude in B-band of 0.0206±0.00I3mag, secondary 
period at 2.05d (or its aliases). 

3.4. LMC164-2 

A x^ test shows that LMG164-2 is significantly variable 
in both ASAS and OGLE datasets. All period searches on 
OGLE data reveal a very clear peak at P = 7.96 d: val¬ 
ues range from 7.9598d for variance methods to 7.9612 d 
for Renson’s method, which is fully consistent with the 
0.0010d Icr error (see Table [3 and left panel of Fig. [3. 
Moreover, this period yields coherent variations with phase 
and, when eliminated from the data, leaves a flat peri¬ 
odogram (Fig. [T] middle and right panels). The EROS-2 
data in blue band (where most data points are available) 
of LMG164-2 have the lowest errors amongst EROS-2 data 
(always<0.045mag). Formally, they do not deviate signifi¬ 
cantly from a constant, but period search methods read¬ 
ily confirm the signal found in OGLE data, with a pe¬ 
riod of 7.959±0.003d and an amplitude in this band of 
0.0089±0.0016 mag. Periodograms and phased variations 
are very similar to those shown in Fig. [3 This signal is thus 
undoubtly real, though it remains buried in ASAS noise. 

3.5. SMC159-2 

SMG159-2 appears only in the OGLE database and the 
OGLE data appear significantly variable, as shown by a x^ 
test. Period searches on OGLE data reveal a very clear peak 
at P = 14.91d : values range from 14.912 d for variance 
methods to 14.915 d for the conditional entropy method, 
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Fig. 5. Period search results for 2dFS936. Left: Comparison of the different period search methods. The presence of a 
signal is marked by a maximum for Fourier and AOV methods, but a minimum for conditional entropy, Renson and 
Lafler & Kinman methods. Middle: Fourier periodogram for the raw (top) and prewhitened (middle, for the best-fit 
Fourier period) data of 2dFS936, along with the spectral window (bottom). Right: Photometric variations of 2dFS936, 
phased using the best-fit ephemeris (see Table [5]). The thick red line corresponds to the best-fit Fourier sinusoid, while 
the symbols for data are as in Fig. [T] 



Frequency 



Fig. 7. Same as Fig. [5] for LMC164-2 (see Tabled for ephemeris). We note in the left panel the presence of peaks at 
frequencies corresponding to multiples of the fundamental period for all methods except Fourier, and in the middle 
panel the best-fit signal at 0.126d“^, its second harmonic at 0.25d“^, and their daily aliases at 0.87d“^ and 0.75d“^, 
respectively. 


a difference which is well within the 0.004 d la error (see 
Table m and Fig. B- This period yields coherent phased 
variations and, when eliminated from the data, leaves a flat 
periodogram (Fig. [5]). Again, this leaves little doubt about 
the reality of the signal. 

Two spectra of SMC159-2 are available, one taken 
in 2013 (on HJD =2 456 640.551) and one in 2014 (on 
HJD =2 456 903.691). Using the best-fit ephemeris derived 
from Fourier algorithm (Table [5|), we found that the 2013 
and 2014 spectra of SMC159-2 were acquired at two very 
different phases, 0.96 and 0.60, respectively. These spec¬ 
tra display many features that are typical of Galactic Of?p 
stars (Fig. [S]). The Hen absorption lines remain essentially 
constant, indicative of a predominantly photospheric origin, 
and they are broader than He i absorption lines, as typically 
observed in Of?p spectra. We note that the apparent spec¬ 
tral type does not change much, the He ii A 4541 /He i A 4471 
ratio remaining consistent with a spectral type of 06. In 
addition, the emission lines appear relatively narrow (e.g. 


FWHM of ~280kms“^ for Ha), i.e. much narrower that 
typical wind Ha emission of 0-stars, as is usually observed 
in Of?p stars. In parallel, many variations typical of Of?p 
stars are also observed. The Balmer lines and Hell A 4686 
emission lines are stronger in 2014, although with one ex¬ 
ception (HD 148937) the differences are generally smaller 
than in the spectra of Galactic Of?p stars. These lines 
appear slightly redshifted at maxi mum, as is seen e.g. in 
GPD-28°2561 (iWade et all. In addition, the Niii 

and Gill emission lines near 4630-4650A are stronger in 
2014, the Gill emissions nearly disappearing in 2013. The 
most extreme variation occurs in the metastable He i A 5876 
line, which changes from an apparent inverse P Gygni pro¬ 
file at minimum, to a much stronger double emission line 
at maximum, the second peak appearing longward of the 
single minimum emission. This could be related to the red- 
shifts observed for other emission lines. If this star is con¬ 
firmed to be magnetic, then the physical origins of these di¬ 
verse behaviours are linked to the geometry of the rotating 
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Fig. 8. Same as Fig. [S]for SMC159-2 (see Table[l]for ephemeris). 
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Fig. 9. Spectra of SMC159-2 recorded in 2013 (^=0.96, thick black line) and 2014 ((/)=0.60, thinner red line) in selected 
spectral windows: clear changes in Balmer and Hei lines can be seen, with stronger emissions in 2014, while photospheric 
Hell absorptions and interstellar lines arising in both the Galaxy and the SMC (e.g. Na D lines near 5900A) remain 
constant. 
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2dFS 936 (P=l/2.48952 d) 



Fig. 6. Photometric changes of 2dFS936, once the 1370 d 
variation is removed. This plot uses the Tq of Table [5] and 
a period of 0.402d (i.e. a frequency of 2.48952d“^). The 
symbols for data are as in Fig. [TJ 
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Fig. 10. Comparison of the peak Ha EW and the dipo¬ 
lar field strength of magnetic O stars (using unpublished 
spectropolarimetry obtained by GAW for NGC1624-2). The 
horizontal red line indicate the peak Ha EW of SMC159-2. 


magnetospheric structure. A full spectroscopic monitoring 
will thus be needed to d erive the stel l ar pro perties, as was 
done e.g. for HD 191612 (jWade et al.Ll201lTl . 

In the Gala xy. Of?p stars bright e n when the ir emis - 
sion increases (IWalborn et al.l . 12004 iBarannik^ l2007ll . 
The phase of the 2013 spectrum of SMC159-2 ((/)=0.96) 
corresponds to a brightness minimum, while the phase of 
the 2014 spectrum ((()=0.60) indicates that it was taken 
close to brightness maximum. Furthermore, a direct com¬ 


parison is possible as the last photometric data were ac¬ 
quired close in time to the 2013 spectrum: the Tmagnitude 
was 15.544±0.003 on i7JZl=2 456 639.571, i.e. the star was 
indeed faint at the time. There is thus a direct correlation 
between the strength of the emission lines and the bright¬ 
ness of the star: SMC159-2 behaves in a manner similar to 
the Galactic Of?p stars. 

Equivalent widths (EWs) of Ha and Hell A 4686 were 
evaluated by integrating the spectra in the 6550-6600 A 
and 4680-4697 A intervals, respectively, as well as by fit¬ 
ting gaussians to the emission peaks. Derived EWs amount 
to 13.7±0.5 A and 4.3±0.2 A in 2013 and to 19.2±0.2 A and 
5.0±0.2A in 2014 for Ha and Hell A 4686, respectively. It 
must be noted that there is no nebulosity near SMC159-2, 
hence the recorded Balmer and He i lines are uncontami¬ 
nated: the emissions and their variability are truly linked 
to the star. We also note that the emission in Ha is very 
strong, suggesting a large quantity of emitting material. 
Amongst magnetic O stars, the peak EW of Ha emission 
gener ally amounts to a few Angstroms (e.g. lHowarth et ^ 
I 2 OO 7 II . except for the exceptional case of the highly magne¬ 
tized NGG1624-2 which presents a peak EW of about 27 A 
(IWade et al.l . I 2 OI 2 DI . Figure [TUI shows the peak EW of the 
Ha line for 0^ OriG and the five Galactic Of?p stars as a 
function of their surface magnetic field strengths B^. This 
figure demonstrates the expected increase in emission with 
the field strength, and hence with the magnetospheric vol¬ 
ume. For slow rotators such as Of?p stars, this magneto¬ 
spheric volume is determined by the Alfven radius: it is thus 
a function of both the magnetic field strength and wind mo¬ 
mentum, but the latter was assumed to be approximately 
uniform for all stars in Fig.jTU] A forthcoming paper (Wade 
et al, in prep) will examine in more detail the theoretical 
expectations for the MG Of?p stars. Based on this figure 
and the outstanding strength of its Ha emission, we can 
predict that the magnetic field of SMG159-2 is very strong, 
second only to that of NGGI624-2 and likely in the range 
5-15 kG. 


4. Conclusions 

The Of?p category comprises magnetic massive stars, 
where an oblique magnetic field provokes rotationally- 
modulated photometric and spectral variations in the opti¬ 
cal, UV and X-ray domains. Up to now, such objects were 
only studied in detail in the Galaxy, but five candidate ex- 
tragalactic Of?p stars have been proposed. No spectropo¬ 
larimetry of these objects is yet available, and further spec¬ 
troscopic data, though somewhat scarce, will be soon pre¬ 
sented by Walborn et al. (in preparation). 

In this paper, we analysed the photometry of these 
five objects. Three of them behave exactly like magnetic 
Of?p stars in the Galaxy: the photometric amplitudes, the 
roughly sinusoidal form of the variations, and the strict and 
stable periodicity with periods from 8d to 1370d are similar 
to what is found in Galactic cases. These stars thus have a 
high probability of being magnetic, though it is difficult to 
derive the system’s geometry and magnetic field strength 
with the photometric data alone. Moreover, simultaneous 
spectral variations observed in SMGI59-2 are also remi¬ 
niscent of the behaviour of Galactic Of?p stars. It can be 
noted that the Ha line of SMGI59-2 is remarkably strong. 
This leads us to predict that its magnetosphere is very large 
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and that its magnetic field is especially strong, likely only 
second to that of the outstanding star NGC 1624-2. 

For the remaining two Of?p candidates, significant pho¬ 
tometric variability is observed, but strict periodicity can¬ 
not be ascertained. Their unexplained photometric varia¬ 
tions may well be a consequence of multiple variability ori¬ 
gins (e.g. rotation in combination with binarity or pulsa¬ 
tion). Therefore, while these stars demonstrate the spectral 
peculiarities expected of hot magnetic Of?p stars, we are 
not yet able to confirm photometric variability consistent 
with stable, long-term rotational modulation. 

This study paves the way for further observations, open¬ 
ing the door to the study of magnetic field and stellar wind 
interactions in extragalactic environment at lower metallic- 
ities. 
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